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ABSTRACT 

We present the discovery of a Type la supernova (SN) at redshift z = 1.914 from the CANDELS 
multi-cycle treasury program on the Hubble Space Telescope (HST). This SN was discovered in the 
infrared using the Wide-Field Camera 3, and it is the highest-redshift Type la SN yet observed. We 
classify this object as a SNIa by comparing its light curve and spectrum with those of a large sample 
of Type la and core-collapse SNe. Its apparent magnitude is consistent with that expected from the 
ACDM concordance cosmology. We discuss the use of spectral evidence for classification of z > 1.5 
SNe la using HST grism simulations, finding that spectral data alone can frequently rule out SNe II, 
but distinguishing between SNe la and SNelb/c can require prohibitively long exposures. In such 
cases, a quantitative analysis of the light curve may be necessary for classification. Our photometric 
and spectroscopic classification methods can aid the determination of SN rates and cosmological 
parameters from the full high-redshift CANDELS SN sample. 



1. INTRODUCTION 

Over the past decade, measurements of Type la super- 
novae (SNe) at redshift z > 1 have extended the observed 
popula tion to a time when the universe was matter dom- 
inated (iRiess et al.l|2001 | |2004 [ J2007| |Suzuki et al.||2012 
Rodney et al.|2012[|Rubin et al.|2U13 1. At these lookback 
times ot > 7 Cyr, the predicted effects of dark energy are 
small, while the typical conditions under which SNe form 
are increasingly different from local environments. 

These characteristics may allow observations at high 
redshift to constrain an evolutionary change in SNIa 
brightness independent of our understanding of dark en- 
ergy. This type of systematic shift in magnitude could 
be ca used by changing metal licity or progenitor masses 
(e.g., Dommguez et al. 20011. Such an effect could be 
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present at a lower level in intermediate-redshift SN sam- 
ples (0.2 < z < 1.0), and therefore be a source of uncer- 
tainty in the determination of the dark energy equatio n- 



of-state parameter (w = P/(pc : ); Riess & Livio 



2006) 



Observations of high-redshift SNe la could also place 
constraints on the binary companions of SN progeni- 
tors. The two most likely SNIa progenitor models are 
the single-degenerate scenario, where a white dwarf ac- 
cretes matter from a main-sequence or giant companion, 
and the double-degenerate scenario, where SNe occur 
through the merging of two carbon-oxygen (C-O) white 
dwarfs. A substantial difference between these mecha- 
nisms, however, is the typical time interval from progen- 
itor formation to explosion; progenitors would likely take 
> 10 9 yr to reach the Chandrasekhar limit by mass trans- 
fer from a nondegenerate companion, but would more of- 
ten take less time in a system of two C-0 whi te dwarfs 
(for a rece nt review of SNIa progenitors, see Wang & 
Han 2012). The distribution of times between forma- 



tion and explosion, known as the delay-time distribution 
(DTD), can therefore be used to set constraints on SN 
progenitor models. Observations of SN rates measure the 
convolution of the DTD with the cosmic star-formation 
history, and h igh-redshift rates are the most sensiti ve to 
delay times ( |Strolgeret al.|2010| [Graur et al.|2011[ ). 

Due to the high sensitivity and angular resolution of 
the Hubble Space Telescope (HST), its Advanced Cam- 
era for Surveys (ACS) has been the most effective instru- 
ment for observing and monitoring SNe out toz« 1.5. 
To find SNe at higher redshifts in the rest-frame optical, 
where they are brightest and we understand them best, 
searching in the near-infrared (IR) with the recently in- 
stalled Wide-Field Camera 3 (WFC3) allows SN surveys 
to reach unprecedented depths not accessible from the 
ground (F160W limiting Vega magnitude ~ 25.5, equal 
to the peak observed brightness of a typical SNIa at 
z « 2.5). 

In this paper we present observations of a SNIa at 



WFC3 F160W + F125W Discovery Image 
2010 Dec 30.7 



5N«: 02:17:46 
6: -05:15:23 





1 arcsec 



Grism dispersion directions 



Fig. 1.— The WFC3 F160W + F125W discovery and difference 
images (using a late-time, SN-free template) for SN UDSlOWil. 
The SN is located ~ 0".l from the center of the host galaxy (2 
ACS pixels). The contours plotted on the difference image of the 
SN (upper left) show the regions containing 68% and 95% of the 
host galaxy light. The center of the nearest neighboring galaxy, 
which causes minimal lensing of the SN (see §4.1), is located ~ 1".5 
away. 

Z = 1.91 (SN UDSlOWil), the highest-redshift SNIa dis- 
covered to date. It was found in the Cosmic Assembly 
Near-infrared Deep Extragala ctic Legacy Survey (CAN- 
DELS, PI: Faber fc Ferguson; IGrogin et al.|2011[ |Koeke-| 
moer et al.||20TT] ). The full CANDELS SN sample is 
designed to measure SN rates and to study SN system- 
atics at redshifts greater than 1. 5. Similar to SN Primo, a 
1.55 WFC3-disc overed SN (|Rodney et al.|2012[|FrecT 



eriksen et al. 2012), UDSlOWil also has spectroscopic 
evidence tor classification. We present the discovery of 
SN UDSlOWil in §2. Section 3 discusses its classifica- 
tion from photometry and HST grism spectroscopy. In 
§4 we estimate the brightness correction due to gravi- 
tational lensing and fit the light curve. We discuss our 
results and the HST spectral classification in §5, and our 
conclusions are given in §6. 

2. DISCOVERY 

SN UDSlOWil was discovered in the second epoch of 
CANDELS obse rvations of the UKIDSS Ultra-Deep Sur 



vey f ield (UDS; Lawrence et al. 2007 Cirasuolo et al 



20071 on 2010 December 30, after subtracting the images 



obtained in the first epoch (2010 November 11). It was 
detected at high significance in both F160W and F125W 
difference images, while a flux decrement was seen at the 
same location in the ACS filter F814W difference image 
(detected at ~ 2.5c). The SN searching is performed by 
eye in the difference images, and in this case we could 
only subtract the first epoch of UDS observations (50 
days before) from the second epoch, as no earlier WFC3 
data were available. The F814W flux decrement suggests 
that pre-maximum SN light was present in the first epoch 
of UDS observations. Thus, the SN was brighter in the 
pre-maximum, shorter-wavelength ACS imaging. 



The WFC3 (F125W + F160W) discovery-epoch image 
of SN UDSlOWil is shown in Figure [l] using a late-time 
(2011 December), SN-free template for the difference 
imaging. The J2000 SN coordinates are a = 02 h 17 m 46 s , 
5 = -05°15'23". It was - 0".l from the center of its host 
galaxy (~ 2 ACS pixels, ~ 0.9 kpc in distance), making- 
it highly probable that this galaxy was the host and un- 
likely that the object was an active galactic nucleus. 

At the time of discovery, we determined the photomet- 
ric redshift of the host galaxy to be > 1.5, although this 
was measured before SN-free WFC3 host-galaxy images 
were available. At this redshift, the SN colors (F160W — 
F814W 3cr upper limit, and F125W - F160W) were con- 
sistent with those expected for a SNIa at z > 1.5 and 
inconsistent with a core-collapse (CC) SN, so we trig- 
gered follow-up observations with the X-shooter spectro- 
graph on the ESO Very Large Telescope (VLT) to obtain 
a spectroscopic redshift of the hostj^j 

Moreover, we monitored the SN with the HST SN 
Multi-Cycle Treasury follow-up program (GO-12099; PI: 
Riess). We imaged the SN with HST (20 orbits, to obtain 
the light curve as well as SN-frec template observations) 
and we obtained G141 grism spectroscopy (15 orbits, for 
resolution R sa 130). 

To measure the photometry of the SN, we subtracted 
the late-time template images from the UDS/SN follow- 
up observations. We measured the flux within a fixed 
aperture of 3-pixel radius and estimated errors in the flux 
from the sky noise of the nearby background-subtracted 
image. Details of the HST observations are listed in Ta- 
ble [ll and the grism spec trum is discussed along with the 
SN classification in 93.21 



2.1. Redshift 

We remeasured the spectral energy distribution (SED) 
of the SN UDSlOWil host galaxy, including photome- 
try from late-time WFC3 and ACS templates as well as 
Subaru, UKIRT, and IRAC data. The Balmer break is 
between the Subaru z band and the WFC3 J band, mak- 
ing the most likely redshift between 1.8 and 2.2 (see the 
lower- left panel of Fig. pjj. Using the X-shooter spec- 
trum, we narrowed this result by identifying [O II] and 
[O III] doublets in the host-galaxy spectrum, giving a 
precise redshift of 1.914. 

The result is also consistent with the HST G141 grism 
spectrum, which contains a clear detection of [O III] 
AA4959, 5007. However, the grism spectrum cannot re- 
solve the doublet, as the spectrum is convolved with 
both the shape of the host galaxy and the point-spread 
function (PSF; combined full width at half-maximum in- 
tensity ~ 116 A) and sampled at a resolution of only 
46.5 A pixel -1 . The VLT spectrum, along with an anal- 
ysis of the late- type host galaxy of SN UDSlOWil, will 
be presented by Frederiksen et al. (in preparation). 

3. CLASSIFICATION 

We classified SN UDSlOWil by analyzing its light curve 
and spectrum, informed by the host redshift. As de- 
tailed below, we first examined the light curve, finding 

15 Based on observations made with ESO telescopes at the 
La Silla Paranal Observatory under program IDs 086.A-0660 and 
089.A-0739. 



TABLE 1 
Photometric Observations 



UT Date 


MJD 


Filter 


Exposure Time 


Flux (counts s x ) 


Vega Mag 


2010 Nov. 


08.8 


55508.1 


F814W 


3517.0 


0.143 ± 0.054 


27.635 ± 0.413 


2010 Nov. 


11.2 


55511.2 


F160W 


1605.8 


0.517 ± 0.074 


25.221 ± 0.156 


2010 Nov. 


11.2 


55511.2 


F125W 


955.9 


0.698 ± 0.096 


25.535 ± 0.149 


2010 Dec. 


28.0* 


55557.4 


F814W 


3817.0 


-0.063 ± 0.041 




2010 Dec. 


30.7* 


55560.7 


F160W 


1705.9 


1.22 ± 0.079 


24.290 ± 0.070 


2010 Dec. 


30.8* 


55560.8 


F125W 


955.9 


1.403 ± 0.102 


24.776 ± 0.079 


2011 Jan. 


12.6 


55573.6 


F160W 


3617.6 


0.901 ± 0.063 


24.616 ± 0.076 


2011 Jan. 


12.8 


55573.8 


F125W 


3617.6 


0.759 ± 0.062 


25.443 ± 0.089 


2011 Jan. 


13.6 


55574.6 


F850LP 


1994.0 


-0.018 ± 0.035 




2011 Jan. 


23.4 


55584.3 


F160W 


3667.6 


0.780 ± 0.061 


24.774 ± 0.085 


2011 Jan. 


23.4 


55584.4 


F125W 


3867.6 


0.535 ± 0.059 


25.823 ± 0.118 


2011 Feb. 


04.2 


55596.1 


F160W 


3767.6 


0.441 ± 0.061 


25.392 ± 0.150 


2011 Feb. 


04.2 


55596.1 


F125W 


3717.6 


0.437 ± 0.062 


26.043 ± 0.154 


2011 Feb. 


16.1 


55608.1 


F160W 


4973.5 


0.309 ± 0.058 


25.779 ± 0.205 


2011 Feb. 


16.3 


55608.2 


F125W 


4973.5 


0.183 ± 0.057 


26.989 ± 0.337 


2011 Jan. 


12.7 


55573.7 


G141 


39088.0 


(grism obs) 




* 

Discovery 


epoch. 













that it is consistent only with a SNIa. In particular, 
the combination of its early-time colors with its rapid 
late-time decline rate does not agree with CC SN mod- 
els. We then used the spectrum to independently rule 
out SNell. While the spectral absorption features alone 
are unable to convincingly distinguish between a SNIa 
and a SNIb/c, SNe II have features t hat are inconsistent 
with the data (see Filippenko 1997 for a review of SN 
spectra). 

3.1. Photometric Classification 

To classify SN UDSlOWil we compared the observed 
UDSlOWil light curve against Monte Carlo simulations 
of Type la and CC SNe at redshift 1.91, generat ed with 
the SuperNo va ANAlysis software (SNANA^J |Kessler 
et al. 2009b). We used a least-squares fit to scale the 
magnitude of the simulated light curves to match our 
data, thus allowing us to examine how our data com- 
pare to the shapes and colors of simulated SNe while 
removing any assumptions on cosmology or intrinsic SN 
luminosity. We then measured the \ 2 statistic for each 
simulated SN compared to our data and converted these 
X 2 values into a Type la SN classification pro bability us- 
ing a simple Bayesian framework, sim i lar to Poznanski 



et al. 



et al. 



Kuznetsova & Connolly (2007), and Sako 
The full Bayesian formalism, along with a 
description of the simulations and our Bayesian priors, is 
presented in the Appendix. 

Our procedure gives us a very high probability that SN 
UDSlOWil is a SNIa. As such, varying our priors on pa- 
rameters such as shape, color, Ay, Rv, or SN rates has 
a very minor effect on the outcome. The reliance on only 
43 CC SN templates constitutes the greatest uncertainty 
in our procedure. Howeve r, using a classific ation proce- 
dure very similar to ours, Sako et al. (2011) found that 
classification using only 8 CC SJN templates still returns 
SNIa classification purities of > 90%. 

We find that the probability of a SNIa was 99.98%, 
with a SNIb/c probability of 2.1 x 10~ 4 (ruled out at 
~3.7ct) and a SN II probability of 1.0 x 10~ 7 (ruled out at 

16 |http : //sdssdp62 . f nal . gov/sdsssn/SNAM-PUBLIC/ . | 



^5. 3c). This indicates that the Type la model dominates 
the probability calculation, and no CC SN models can 
adequately describe the observed photometric data. 

Figure [2] shows the best-fit light curves, along with 
the flux range of simulated SN light curves encompass- 
ing 95% of the Bayesian evidence for each SN type. The 
best-fit light curves for Types la, Ib/c, and II SNe re- 
turn X 2 /v values of 18.6/11, 35.5/11, and 51.1/11, re- 
spectively. Note that these x 2 values are only illustrative 
of the quality of the match for each model. They rep- 
resent the best match from a large but limited number 
of random simulations, so one cannot use these values in 
X 2 goodness-of-fit tests for model rejection. By contrast, 
the final classification probability relies on the weight of 
evidence from all realizations of each model. 

Our best-fit x\ and C values for the Type la model 
were —1.56 and —0.12, respectively. These values are 
fully consistent with t he SA LT2 parameters derived from 
light-curve fitting in £4.2 (x\ = —1.50 ± 0.51 and C = 
— 0.07±0.11). We note that if we increase the errors such 
that the SNIa x 2 l v ~ 1 (accounting for the possibility 
that we underestimated the uncertainties), the Type la 
probability is still as high as 99.84%. Figure [2] shows that 
the nearly 100% probability of classification as a SNIa 
(and the superior best-fit x 2 value) arises because the 
SNIb/c and SNII light-curve fits are unable to match 
the combination of SN UDSlOWil's high signal-to-noise 
ratio (S/N) discovery-epoch colors and its rapid light- 
curve decline rate. 

As a verification of this light-curve classification, we 
used the Photometric Sup erNova IDcntification software 
(PSNID; [Sako et al.|[2fJ08| , finding that it also prefers a 
SNIa with a slightly higher 4.1a confidence. The differ- 
ence in probability is primarily due to our conservative 
CC model uncertainties (see the Appendix), which re- 
duce the x 2 values of CCSNe. Although it only uses 
8 CC SNe, the purity of PSNID classifications has been 
robustly tested using Sloan Digital Sky Survey (SDSS) 
SNe, and it obtained the highest figure o f merit in the 
SN P hotometric Classification Challenge (Kessler et al.| 



2010) 
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Fig. 2. — The best-fit light curves using SNANA-simulated Type la, Type Ib/c, and Type II SNe (top to bottom). Red and green colors 
indicate photometry and simulations for the F160W and F125W filters (approximately rest-frame B and V), respectively. Shaded regions 
include the range of SN fits that encompass 95% of the Bayesian evidence for a given SN type. Type Ib/c and Type II fits have higher x 2 
values (even including model uncertainties), due to their inability to match the combination of the colors nearest to maximum light and the 
rapid rate of decline. The sum of the evidence for a SN la model gives a nearly 100% probability that the SN is of Type la. We included 
F814W and F850LP data in the \ 2 fitting, but have omitted them in this figure for visual clarity. The fluxes shown are transformed from 
Vega magnitudes using a zeropoint of 27.5. 

from the grism spectrum. We omitted these pixels (the 
shaded region in Fig. [3]) when we later fit spectral tem- 
plates to the SN spectrum, as we did not have a SN-frcc 
line measurement to subtract from the observations. We 
then rescaled the aXeSim output spectrum to match the 
F160W magnitude of the host galaxy as measured in the 
last epoch of follow-up imaging after the SN had faded. 
Our simulated host-galaxy spectrum is shown in Figure 
3] (upper left). 

After subtracting the host-galaxy model from the SN 
spectrum contaminated with host light, we u sed the Su- 



3.2. Spectrum 

Spectroscopic confirmation of SNe has pr oven challeng- 
ing at the se redshifts (see the discussions in Rodney et al. 
2012| and [Rubin et al.|[2013| , due to the difficulty of ob- 
tammg high S/JN observations and the paucity of defining 
features in the available window (for UDSlOWil, ~ 1.12— 
1.65 /xm; rest frame ~ 3840-5660 A). In the case of SN 
UDSlOWil, the SN was separated from its host galaxy 
by only ~ 0".l, contaminating the SN spectrum with 
host-galaxy light. We removed the host galaxy from the 
spectrum by subtracting a section of the galaxy free from 
SN light, but the combined noise from the SN and host- 
galaxy spectra made a spectral classification inconclu- 
sive, even with substantial host-galaxy smoothing. 

As an alternative approach that avoids adding addi- 
tional host-galaxy noise to the SN spectrum, we gener- 
ated a noise-free synthetic host spectrum. We fit SEDs, 
using a library of spectral templates, to optical and near- 
IR Subaru, ACS, WFC3, and UKIRT host-galaxy pho- 
tometry following |Dahlen et al. (2010). We then simu- 
lated the observed grism host spectrum with the aXeSim 
software package p^ The aXeSim software convolves the 
SED with the shape of the host galaxy and HST PSF and 
samples the spectrum at the G141 spectral resolution of 
46.5Apixer 1 . 

One would not necessarily expect emission lines to be 
the same strength in the template as in the real galaxy 
due to its differing metallicity, star formation rate, and 
population of massive stars. Therefore, we replaced the 
pixels covering the [O III] line in our simulated host 
galaxy with those covering the prominent [O III] line 



1 



pernova Iden tification (SNID) codq | from Blondin & 



http : //axe . stsci . edu/axesim/ . 



Tonry (2007) to match the UDSlOWil spectrum with 
Type la, Type Ib/c, and Type II SN template spectra 
to determine the best-fit spectrum for each class. For 
SNIa fits, we only used templates within ±3 rest-frame 
days of the age of the SN UDSlOWil spectrum. The age 
is based on the SALT2 fit in §4.2, which gives ~ 12 ± 1 
day after maximum (rest frame). For CC SN fits, we used 
any templates which put the time of maximum between 
the first two epochs of observation (the same as our prior 
in §3.1). When fitting the spectrum to templates, SNID 
removes the continuum using a high-order polynomial fit 
and only matches the spectral features, making the fit 
independent of reddening and brightness. 

SNID returns the rlap parameter, which is meant to 
quantify th e quality of the correlation (see 
Tonry|2007| for details) " 

that rlap values less than 5 are inconclusive. JNote that 
SNID does not apply a broadening symmetric function 



18 http : //www . oamp . f r/people/blondin/sof tware/snid/ 
index.html. 
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Fig. 3. — The spectrum of SN UDSlOWil, from which the host-galaxy light has been subtracted using a spectrum fit to the host's SED. 
The host SED and its best-fit spectral template is shown at the bottom left. The two-dimensional and one-dimensional grism spectrum 
along with the simulated host-galaxy spectrum (the SED convolved with the shape of the host using aXeSim) are given at the upper left. 
On the right we show the best-fit SN templates from SNID for Types la, Ib/c, and II. We illustrate median bins to visually emphasize the 
spectral features, but have used the unbinned data for our analysis. Although the spectrum can be fit reasonably well by both a SNIa or 
a SNIb/c, a SNII does not contain the spectral features seen at ~ 4600 and ~ 5200 A in the rest frame. Data in the [O III] subtraction 
region (the shaded region) were not used in the fit. 



(Tonry & Davis 1979), which should be used when the 
widths of SIN features are comparable to the resolution 
of the spectrum. As the rest-frame G141 resolution is 
~ 16 A pixel -1 and significant SN features have a typ- 
ical width ~ 50 A, this is not a major concern. How- 
ever, the inclusion of this function for grism data could 
improve future SNID classifications, especially those at 
lower redshift. 

The right side of Figure [3] shows the best-fit Type la, 
Ib/c, and II SN templates with rlap values of 4.8, 2.7, 
and 1.9, respectively. We show median bins to emphasize 
the spectral features, but we fit spectra to the unbinned 
data. The data can be fit by 6 other normal SNe la with 
rlap values of at least 4. Five other SN Ib/c fits have an 
rlap of at least 2, and only two other SNII fits have an 
rlap of greater than 1.5. 

The x 2 values for the fits, now including continuum, 
to Type la, Ib/c, and II SNe are 95.9, 97.9, and 104.0 
(respectively) with 106 degrees of freedom. We note that 
both SNe la and SNe Ib/c can provide good fits, although 
the former give a slightly better match. However, all 
SNII templates yield a poor correlation; the rest- frame 
features at ~ 4600 and ~ 5200 A (which are created by 
neighboring Fe II, Fe III, Si II, and S II absorption in 



Filippenko 1997) are not well fit by the spectra 



SNe la; 

of SNell. We note that although not all of the SNID 
SNe II are as featureless as the best-fit spectrum shown 
in Figure [31 they all have difficulty matching the strength 
or location of the spectral features. Even Type II-P tem- 
plates, which typically have stronger features, have a 
maximum rlap of only 1.4. 

Because the best rlap value is less than 5, this SN 
cannot be considered to be spectroscopically confirmed. 
In addition, the "lap" parameter, which describes the 
overlap in wavelength space between the template and 
SN spectrum, is 0.39 for each of the best-fit SNIa and 
CC SN t emplates. This is below t he minimum lap of 0.4 
used by Blondin & Tonry (2007) for spectral confirma- 
tion. However, the SIN UDSlOWil spectrum still favors 
classification as a SN la and its rlap is com parable to that 
of oth er high-redshift SNe la. SN Primo ( Rodney et al.| 
z = 1.55) had an rla p of only 3.7. \Ve also fit the 
spectrum of SN SCP-0401 (|Rubin et al.|2013| z = 1.71), 



finding that it is best matched by featureless SNID Type 
la and Type Ib/c spectra. However, if we require a match 
to at least one spectral feature or a lap value greater than 
0.1, the maximum SNIa rlap is 4.6 (SN 1993ac, +7 days; 
lap of 0.2). SNID templates begin at A rcst = 2500 A, so 



G 



the first ~ 500 A of the SN SCP-0401 spectrum were not 
included in the SNID fit. 

4. ANALYSIS 

Taken together, the photometric evidence suggests 
that UDSlOWil is a SNIa with high confidence. The 
spectroscopic evidence independently favors a Type la 
classification. We now proceed to derive its shape and 
color-corrected magnitude, taking into account the pos- 
sibility that the SN light has been gravitationally lensed 
by foreground structure. 

4.1. Lensing 

Our ability to use SNela as accurate distance indi- 
cators to constrain cosmological parameters requires us 
to determine the impact of foregro und-matter in homo- 
geneities on the flux of the SN (e.g., Jonsson et al.|2006). 
Even without multiple images, gravitational lensing can 
significantly magnify the SN, altering our measurement 
of its distance. SN UDSlOWil is close in projection to 
another galaxy (see Fig. [I]) separated by only 1".54. 
Therefore, it is necessary to estimate the possible mag- 
nification of the SN which could lead to a bias in the 
derived SNIa distance. All other foreground sources are 
greater than 4". 5 away and cause negligible magnifica- 
tion. 

We fit the SEP of the lens galaxy as described by |Wik- 
lind et al. (2008) to characterize its physical properties. 
We u sed a (Jhabrier initial mass function (IMF; Chabricr 
20031) rather than the Salpeter IMF cited by |Wiklind| 
et al. (20081; the former gives a slightly smaller stel- 
lar mass but is a more accepted model. To account for 
photometric uncertainties, we drew Monte Carlo sam- 
ples for the measured photometry of the galaxy and 
used the best-fitting SED models to characterize the 
SED. The SED fit indicates a low-mass galaxy with a 
photometric redshift 0.283 ± 0.080 and a stellar mass 
log(M*/M ) = 7.968 ± 0.222. We used these parame- 
ters to create a plausible mass model of the galaxy and 
estimate the magnification of the SN. 

We modeled the stellar component of the galaxy as an 
expon ential disk u sing parameters measured from GAL- 
FIT (Peng et al. 2002[ ) and the dark-ma tter halo using 



a Na varro-Frenk- White (NFW) profile (Navarro et al 



f J7t. We used the broken power law given by Yang et al. 
T2| to relate the stellar mass M* to the halo mass and 
mass-concentration relation given by jMaccio et al.| 

(2008) when modeling the halo. 

Both the mass-concentration and stcllar-to-halo mass 
relations have significant scatter around the median rela- 
tions. To account for this scatter, we took 10,000 Monte 
Carlo realizations of lensing potentials to calculate the 
expected magnification distribution. We also drew a 
photo-z and stellar mass distribution from the Monte 
Carlo realizations of our SED fits. 

Despite the proximity of the galaxy, its low mass makes 
magnification a minor effect. The median magnification 
from the above analysis is 2.8JT^ '|%, where the lower 
and upper uncertainties represent the 16 th and 84 th per- 
centiles (±1ct), respectively. These models do assume a 
spherical NFW profile, but adding ellipticity to the halo 
does not significantly change our results. This analysis 
shows that the systematic offset due to lensing is much 



smaller than the photometric uncertainties; we applied 
the lensing correction to our derived magnitude, but it 
does not have a significant effect on our distance modu- 
lus. 

4.2. Light-curve Fit 

We fit the light curv e using the SALT2 implementa- 
tion ( |Guy et al.||2010h contained in SNANA (Fig. E). 
Note that although the ACS data provide a valuable 
color constraint for classification purposes, we have omit- 
ted them from our cosmological analysis. At z = 1.914, 
the ACS bands sample rest-frame wavelengt hs of 2400- 
3300 A, where SNela are more hetero geneous ( Ellis et al 
2008|) and may evolve with redshift (IFoley et al.||2U12) 



Furthermore, the rest-fram e UV has been prob lematic 



for SN la light-curve fitters ( Kessler et al. [2009a ) . Given 
these concerns, we discarded the ACS data tor our light- 
curve fitting in order to avoid introducing a bias in the 
derived distance. The ACS observations provided only a 
single measurement with positive flux (F814W at MJD 
= 55801.1), so this does not exclude a large fraction of 
useful data. 

The light-curve parameters for SN UDSlOWil are typ- 
ical of SNela; we derive values of x\ = —1.50 ± 0.51 and 
C = —0.071 ± 0.11, which are cons istent with the pa- 



rameters described by 
a c =0.13, xi = -0.13, 
a peak magnitude m* B 



-- 0.04, 
TA also gives 
26.20±0.11. We then converted 



our SALT2 values into SiF TO values dCo nley et al 

SN 



using the relations of |Guy et al. (20 



2008J), 
in order to use 



the shape and color c onstants "from sNLS (a = 1.367, 
/3 = 3.179; |Sullivan et al.|20lT] >. We derived a light-curve 
shape and color-corrected magnitude (m cor r) using 



+ a x (s - 1) - x C, 



(1) 



where m cmr is equal to the distance modulus plus the 
SN absolute magnitude, M. Here m* B is the peak SN 
magnitude, s is the SiFTO stretch parameter, and C de- 
scribes the color; also, m* B includes the lensing correction 

of 0.030t^i3 ma S- 

This analysis gives a corrected magnitude of 26.15 ± 
0.39. We compare to the corrected magnitude for 
ACD M, m\cr)M, by using the cosmological parameters 
from |Sullivan et at] p01l| ) (fi A = 0.73, VL M = 0.27, 
71.6kms _i Mpc -1 ) and a least-squares 



1, H, 



o 



fit to the |Conleyetal.| ( |2011[ ) SNe. We added an offset of 
0.27 mag to the value of m B for the Conley et al. (2011 1 
SNe in order to match the normalization of the SAL'1'2 
fitter contained in SNANA, finding tukcdm = 26.46 mag 
(including the offset, this gives an absolute SN magni- 
tude M = -19.39). SN UDSlOWil is less than la from 
ACDM. 



We also fit t he light curve with MLCS (Jha, Riess, & 
Kirshner|2007 l, after using SDSS SNe to determine, and 
correct for, the TO corr offset between MLCS and SALT2 
fits. MLCS gives the same corrected magnitude with a 
somewhat smaller uncertainty, TO corr = 26.15 ±0.27 mag 
(with A = 0.30 ± 0.18 and A v = 0.01 ± 0.05 mag). This 
value is slightly brighter than expected from ACDM, but 
consistent at 1.15c. We verified that the MLCS and 
SALT2 param eters are cons i stent w ith each other using 



relations from Kessler et al. (2009a) 
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Fig. 4. — On the left, the SALT2 light-curve fit to SN UDSlOWil. SALT2 gives a normal set of light-curve parameters along with a 
corrected magnitude of 26.15 ±0.39 (consistent with ACDM) and a reduced \ 2 of 1.5. On the right, we place SN UDSlOWil on the Hubble 
residual diagram using cosmologica l parameters from [Sullivan et al.| || 201lj| and Ho = 71.6k ms —1 Mpc - 1 . For comparison, we also show 
the compilation of ~ 500 SNe from [Conley et al.| | |201l| |. .Lastly, using |bullivan et al. ] ( |2011[ | values for a. and 0, we include the othe r two 
opic evidence for classification discovered at a redshift greater than 1.5 l |Rodney et al.|2012||Rubin et al.|2013fl . The 
1 ''"' ' ' i • i '"- ick collaboration (M\ = 0.685, 



SNe la with spectroscopic < 



dotted line indicates the differenc e in m\ar>M when using the recent cosmological parameters from the 
Q M = 0.315; |Planck Collaboration et al.|2013[ |. 
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At z > 1.5, deriving a spectral classification with HST 
requires a large number of orbits to obtain a high S/N. In 
addition, the HST IR grisms cover a relatively small rest- 
frame wavelength range. In the case of SN UDSlOWil, 
the G141 grism wavelength range (~ 1.12-1.65 /im; rest 
frame ~ 3840-5660 A) does not include either the Si II 
absorption at ~ 6150 A or the Ca II H&K trough at ~ 
3750 A, which are some of the deepest SNIa features (the 
features have equivalent widths ~ 100 A). This means 
that spectral classifications of high-redshift SNe using 
SNID will have lower-significance correlations with SN la 
template spectra, and thus often yield rlap values much 



less t han the suggested minimum of 5 (Blondin & Tonry 



1.9 

redshift 



Fig. 5. — The number of orbits necessary to rule out the pos- 
sibility of a SN II or SN Ib/c 90% of the time when observing a 
SN la. Using aXeSim, we simulated a variety of exposure times in 
the redshift range 1.5-2.3. We found that ~5— 10 orbits can rule 
out a SN II, but that the number of orbits required to rule out a 
SNIb/c possibility is significantly greater. The number of orbits 
to rule out a SNIb/c is much lower in the region where Si II and 
Ca II H&K are completely visible. These results indicate the need 
for photometric evidence in SNIa classification at high redshift, 
although the HST grism can also be valuable in determining SN 
redshifts. 

5. DISCUSSION 

The observations of SN UDSlOWil presented here 
demonstrate that the HST WFC3 now allows the cos- 
mological study of SNe la at higher redshifts than ever 
before. The analysis presented above is enabled by the 
photometric classification methods we employ. However, 
SN science, especially cosmology, has in the past relied 
heavily on spectroscopic evidence for classification. 



20071. 

Host-galaxy contamination can be a significant source 
of noise in high-z SN grism spectra. In this work we 
have used aXeSim to remove host-galaxy light from the 
SN+host spectrum. However, even in a situation where 
a SN is well separated from its host galaxy, spectral evi- 
dence alone may not be enough to unequivocally classify 
the SN as Type la. Figure [5] shows the number of HST 
orbits with the G141 IR grism that are needed for SNID 
to correctly distinguish a SN la from a SN Ib/c or a SN II 
90% of the time, in the redshift range 1.5-2.3. For this 
figure, we used aXeSim to simulate 100 SNIa observa- 
tions per unit redshift with 2.6 ks per orbit, from to 100 
orbits (2-orbit intervals). The background flux was set 
to 95 e~ s _1 , and the simulated SN magnitude was fixed 
at the peak magnitude of SN Primo (an F160W Vega 
magnitude of 23.98; Primo was observed closer to maxi- 
mum than UDSlOWil). Each simulated grism spectrum 
was then processed with SNID, and the classification was 
deemed correct if the best rlap for a SN la template was 
larger than the best SN Ib/c or SN II rlap (we allow rlap 
<5). 

Figure [5j shows that, similar to the case of SN 
UDSlOWih ruling out a SNII possibility requires only 



~5-10 orbits (1.5 < z < 2.3). For a large program like 
CANDELS, with 200 follow-up orbits and a desired sam- 
ple of ~10 SNe, this is a feasible number. Ruling out a 
SN Ib/c, however, can require up to 25 orbits, becoming 
most costly at those redshifts where Si II and Ca II H&K 
are not completely visible. At very high redshift, such as 
z « 2.2-2.3, Figure [5] shows that the G141 exposure time 
required to distinguish a SN la from a SN Ib/c begins to 
drop. The value of spectroscopic confirmation for such 
high-redshift SNe may warrant the necessary investment 
of orbits, especially if additional high-value targets can 
be simultaneously observed within the grism field of view. 
We note that simulating CC SN observations shows that 
SNID can occasionally misclassify CC SNe as SNe la, an 
effect we have not taken into account in this analysis. 

The G102 grism can also be useful for picking out 
features such as Si II and Ca II H&K in the redshift 
ranges where the G141 grism does not contain them. 
Unfortunately, the consequence of its more limited wave- 
length range (~ 0.8-1.15 /zm; rest frame ~ 2750-3950 A 
at z = 1.91) is that a SNIb/c template is more likely to 
match a SN la G102 spectrum well. 

HST grism spectroscopy can be good at determining 
SN redshifts or host-galaxy properties. However, for reli- 
able SN classification, photometric evidence is often im- 
portant. In the case of SN U DSlOWil we improve d upon 



and its light-curve fit (with SALT2) is consistent with 

T 



ACDM. An alternative fit with MLCS (Jha et al. 20071 



the photometric methods of Rodney et al 



(20121 by in- 



troducing a quantitative Bayesian method that returns 
probabilities for each SN type. With both photometric 
and spectroscopic methods, we can be confident in our 
classification and subsequent analysis. 

6. CONCLUSIONS 

At a redshift of 1.914, SN UDSlOWil is the most dis- 
tant SN la yet known. Classification of this SN rests on 
photometry and grism spectroscopy, which rules out the 
possibility of a CC SN. The spectral evidence alone disfa- 
vors the possibility of a SN II, while supporting a SN la or 
SNIb/c hypothesis. The combined SN colors and rapid 
decline rate are inconsistent with a CC SN and in good 
agreement with a SN la model. 

We find that SN UDSlOWil is not significantly lensed, 



is slightly brighter than ACDM, but consistent at 1.15<7. 
When the full analysis of the CANDELS SNe is com- 
plete and combined with the data from the Cluster Lens- 
ing and Supernova survey with Hubble (PI: Postman; 
Postman etal.|2012[ ), we expect that SN UDSlOWil will 
be one of a sample of ~ 10 SNe la above a redshift of 1.5 
to be found by these programs. This SN is an example 
of an object in a new area of SN cosmology, one which 
has only begun to be explored in the last few years with 
the advent of WFC3 on HST and one with unique classi- 
fication challenges. However, with the full sample of SNe 
at redshift greater than 1.5, new limits on the evolution 
of dark energy, the DTD, and the evolution of the SNIa 
population will become possible. 
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APPENDIX 
PHOTOMETRIC CLASSIFICATION METHOD 



We began our classification procedure by using SNANA ( Kessler et al.|[2009b[ ) to gener ate a Monte Car lo simulation 
of 30,000 SNe at redshift 1.91. 10,000 simulated SNe la were based on the SALT2 model ( |Guy et al.|2010[ ), with values 
of the shape parameter x\ drawn uniformly in the range —3 to 3 a nd the color paramet er (J from —0.4 to 0.6. These 
ranges cover the observed distribution of SALT2 pa rameter values ( Kessler et al.|2009a I , and the C term accounts for 



both intrinsic SN color and host-galaxy extinction ( |Guy et al.||2007 ). The remaining 20,000 simulated SNe were split 
evenly between the two principal CC SN classes, with light curves based on 16 Type Ib/c and 27 Type II SN templates 
that comprise the SNANA non-SNIa library (including subtypes lb, Ic, II-P, Iln, and II-L). Host-galaxy reddening 
was applied to each simulated CC SN using Ry — 3.1, with a random draw of Ay in the range to 7 mag using the 
Cardelli et al. ( 1989| reddening law. For both the SNIa and the CCSN simulations we chose random values for the 
date of the light-curve peak, using a range spanning the first to the second epoch of UDSlOWil observations. 

To compare each of the 30,000 synthetic light curves to the N = 15 photometric observations of SN UDSlOWil, we 
computed the \ 2 statistic given by 



* 2 = E 



V*obSji -"• - x - ™sim,i) 



(Al) 



obs,, 



where -F bs,j and <7 bs,i are the fluxes and uncertainties for each observation. Here .Fsim.j and er s im,i are the fluxes and 
uncertainties (respectively) for a single simulated SN on each observation date. The variable A is a scaling parameter, 



described below. For SNela, most of their intrinsic variability can be described by the SALT 2 model's shape a nd color 



parameters. Additional variability causes scatter about the Hubble diagram, and is given by Guy et al. ( 2010| as 8.7% 
in distance modulus. We treat this variability as approximately equal to the model uncertainty, which m hux space 



translates to er s i m i 



Simla, i 



0.08,4 x F si 



CC SNe have greater heterogeneity, such that our relatively small set of discrete templates cannot describe the 
entire population. By setting a nonzero er s imcc,ij our limited CC SN temp late library can more accurately represent 
this diverse class. Considering a similar problem, Rodney & Tonry (2009) estimated <7 s imCC,i by measuring the flux 
difference between all possible pairwise comparis ons ot templates of the s ame subclass and taking the median. We 
classify SN UDSfOWil using more templates than Rodney & Tonry (20091, such that the CCSN population is better 
sampled and less uncertainty is present. However, we adopt their value ot cr s ; m cc,i = 0.15-4 x F s i mj i as a conservative 
estimate. 

We next chose the optimal distance or absolute magnitude of every simulated SN, therefore removing the assumptions 
on cosmology and SN luminosity functions that are built into the SNANA simulations. Here we have multiplied -Fsim.i 
by A, a free parameter that introduces a coherent flux scaling across all bands. We find a separate value for A with each 
of the 30,000 comparisons, using x 2 minimization to match the simulated magnitudes to the data with the equation 
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We then converted the measured \ 2 value s into a Type la SN classificati on pr obability us i ng a simple Bayesian 



framework, similar to Poznanski et al 
that the data (D) match a simulated fcW 
types, and time of maximum light) is given by 



( 2007[ ), |Kuznetsova fc Connolly| ( |2007| , and |Sako et 
N of type 1 j with parameters H (shape xi, color C f< 



al. (2011). The likelihood 



ToFSNTa~6r A v for CC SN 



p(D\d,Tj) 
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where x 2 is given in Equation 1. Multiplying by prior probability distributions for each of the model parameters then 
gives us the posterior probability for each point in parameter space, p{9\Tj)p{D\9 , Tj). As we are interested in model 
selection, not parameter estimation, we can marginalize over all of the nuisance parameters 6. Approximating the 
marginalization integral with a discrete sum, the probability of SN type Tj given the model is 



JV Blm (T 3 ) 



p{D\T j )= y, pWTMmT^Sd. 



(A4) 
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For SNIa parameters X\ and C, we applied Gaussian priors based on the values given by Kessler et al. (2009a) 
(C = 0.04, erg = 0.13, xi = — 0.1 3, a Xl = 1.24). For t he CC SNe pa r amete r Ay, we used the Monte Carlo recipe 
provided by Riello & Patat (2005) and implemented by Dahlen et al. (20121 for a random galaxy orientation. The 
distribution is peaked at Ay = mag and falls off quickly such that Ay ■> 3 mag is very unlikely. We used a flat prior 
for the time of peak brightness. 

Note that for computational efficiency we have used SNANA to sample the multi-dimensional model parameter space 
using a Monte Carlo simulation with uniform sampling distributions (instead of the more typical approach, using a 
grid of parameter values). Thus, we must approximate 56 - the vector of step sizes along each dimension of parameter 
space - using the range over which each parameter is sampled: 
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where Adk is each range, Nj is the number of simulated SNe in the class (we used 10,000), and the product is over 
Ng, the number of parameters 9 for the model: 3 for SNela (xi, c, t p k) and 2 for CC SNe (Ay, t p k). 
Lastly, we multiplied each model by a SN rate prior P(Tj). This prior is the fraction of SNe at re dshift 1.91 that we 



expect to be a given type. We began by adopting the rate measurement from |Dahlen et al. ( 2008 ) for SNe la and the 



of|Ho pkin 

2oo3rrr 



lo cal rates from Li et al. (|2011 1 for CC SNe. W e scaled the CCSN rates according to the cosmic star-formation history 



ins & Beacom 1 200t> ) using the form of Cole et al. (2001) and a modified Salpeter IMF (Baldry & Glazebrook 



2003p. The normalized rates showed, as an estimate, that one could expect only ~ 2% of SNe at this redshift to be ot 



1'ype la. The SN UDSlOWil host galaxy SED (§3.2) is consistent with a starburst galaxy, so it is possible that these 
ave rage rates ov erestimate the SNIa rate in this galaxy. In addition, SN rates are very uncertain at this redshift, and 
the Dahlen et al.| (|2008[) rates at this redshift are consistent with SNe la; however, we note that lowering this rates 



10 



prior by an order of magnitude still returns a classification probability greater than 99%. 
independent of the uncertainty in SNIa rates. 
Applying Bayes' theorem gives the final probability that SN UDSlOWil is of Type la: 



p(Ia\D) = 



p(D\Ia)p(Ia) 

Ep(d\tMt 3 [ 



Thus the result is largely 



(A6) 



where the summation is over Type la, Ib/c, and II SN models (Tj). In the case of SN UDSlOWil, our likelihood 
function is sufficiently narrow that the priors have only a minor effect. Thus, we found that allowing Ry, SALT2 
parameters a and /?, or parameter ranges to vary does not substantially alter the high probability that this SN is of 
Type la. 
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